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C/3 ' Abstract. We study the response of single nanosized spherical colloids in electrolyte 

solution to an alternating electric field (AC field) by computer simulations. We use a 
coarse-grained mesoscopic simulation approach that accounts in full for hydrodynamic 
and electrostatic interactions as well as for thermal fluctuations. The solvent is modeled 
'^ I as a fluid of single Dissipative Particle Dynamics (DPD) beads, and the colloidal 

C ' particle is modeled as a rigid body made of DPD beads. We compute the mobility 

fj . and the polarizability of a single colloid and investigate systematically the effect of 

amplitude and frequency of the AC-fields. Even though the thickness of the Debye 
layer is not "thin" compared to the radius of the colloid, and the thermal fluctuations 

^ , are significant, the results are in good agreement with the theoretical prediction of the 

CO ' Maxwell-Wagner-O'Konski theory, especially for uncharged colloids. 
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1. Introduction 

Alternating electric fields (AC fields) can be used to manipulate single colloid or colloidal 
suspensions. This provides an attractive tool for controlling the position of individual 
colloid or the structure of colloidal assemblies. Since the colloids respond to external 
fields on relatively short time scales and in an often fully reversible way, the approach 
provides high flexibility in designing optimal AC pulses for a given task. 

One particularly prominent example is dielectrophoresis. In dielectrophoresis, one 
exploits the fact that external electric fields polarize the colloids, which are then in turn 
driven along the direction of the field gradient. The time-averaged force acting on a 
particle in the AC- field 'Eexp{iut) is 

Fdep ~ Re{a(u;)}V|E|^ (1) 

where a(uj) is the complex polarizability of the particle. The dielectrophoretic effect 
can be used for trapping colloidal suspensions in "electric bottles" [H [2], but also for 
selectively trapping single colloid and separating colloids or viruses by size and shape 
[3111]. It has thus many potential applications in nano- and biotechnology. 

The crucial quantity in dielectrophoresis is the polarizability a{u). Assuming that 
the external field has a form Eexp{iut), the effective dipole moment of a spherical 
particle can be written as 

p = a(w)E = 4:7ieraK{uj)R^ E, (2) 

where e^ is the permittivity of the surrounding medium, K{u) is the Clausius-Mosotti 
factor and R is the radius of the particle. For uncharged particles, the polarizability is 
governed by the permittivity and conductivities of the particle and the medium. For 
charged particles in salt solution, the electrical double layer surrounding the particle 
contributes significantly to the polarizability, and the dielectrophoretic response results 
from an intricate interplay of dipole reorientation, ion mobility, electrostatics, and 
hydrodynamics. The dielectric response of colloids to AC fields has attracted abiding 
interest for decades, but it is still not fully understood, especially in cases where the 
colloids have nanometer size and the electric double layer is thick, since most theoretical 
studies neglect thermal fiuctuations and assume that the electric double layer is thin 
compared to the colloidal radius. 

Several processes, on time scales spanning several orders of magnitude, contribute 
to the dielectric response of charged colloids in salt solutions. The fastest event following 
the sudden apphcation of a constant electric field is the interface polarization induced 
by the local dielectric contrast between water and the colloid material. This happens on 
the GHz scale (the relaxation time of water is 17 GHz), which is instantaneous for the 
purpose of the present discussion. Likewise, the time for hydrodynamic flow to build 
up is also almost instantaneous [5]. On the MHz scale (1-10 MHz), ionic migration in 
the electrical double layer sets in and dominates over the molecular dipole reorientation. 
In this regime, the main contribution to the induced dipole moments stems from the 
conductivity mismatch between the particle and the solvent due to the presence of free 
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counterions in the solvent (Maxwell- Wagner relaxation). Finally, on long time scales up 
to lOOHz, a salt concentration gradient builds up along the colloid and the thickness 
of the double layer varies accordingly, leading to an additional source of polarization 
(a-polarization). 

Various approaches have been proposed to explain the dielectric response of colloidal 
particles. In the high-frequency regime, the Maxwell- Wagner mechanism of dielectric 
dispersion [6l |7] has been widely discussed, which only depends on the bulk properties 
of the solution and the colloid. O'Konski later introduced a surface conductivity term 
to account for the contribution of the electric double layer [8]. The Maxwell- Wagner- 
O'Konski (MWO) theory has also been extended to ellipsoidal particles [9]. Theories 
for the low-frequency region have been developed in the Ukraine school [10] based on 
the standard electrokinetic model [H]. Most studies rely on the assumption that the 
electrical double layer is much thinner than the radius of the particle. For situations 
that involve thick electrical double layers and the whole frequency spectrum, numerical 
methods to solve the electrokinetic equations have been implemented p^ [T3l ITi] . 
In the standard electrokinetic model, ions are treated as point particles, the size 
effect is neglected, and a mean field approach is taken, i.e., thermal fluctuations are 
neglected. Moreover, linear response theory is employed in the numerical calculations, 
but experiments have suggested that the dielectric response becomes nonlinear already 
for moderate field strengths |15j . 

Molecular simulations can shed light on dielectrophoretic phenomena in a well- 
defined model system, but are still scarce [Ml [17]. Such studies are numerically 
challenging, because two different types of long-range interactions are involved: the 
electrostatic and the hydrodynamic interactions. In recent years, a number of coarse- 
grained simulation approaches have been developed to address this class of problem. The 
idea is to couple explicit charges (ions, colloids) with a mesoscopic model for Navier- 
Stokes fluids. On the side of the fluid model, different approaches exist in the literature 
that range from Lattice Boltzmann (LB) methods [TH [19], Multi-Particle Collision 
Dynamics (MPCD) [201 [21], to Dissipative Particle Dynamics (DPD) [221 [231 [21]. In 
this paper, we use the DPD method. DPD is a coarse-grained simulation method which 
is momentum-conserving and Galilean invariant. Since it is a particle-based method, 
charges can be introduced quite naturally. A recent comparative study [25] indicated 
that the Coulomb interaction is the most time-consuming part in the simulation. The 
computational costs of different methods for modeling the fluid thus becomes comparable 
at intermediate or high salt concentrations. 

In the present paper, we use DPD simulation to study the dielectric response of a 
single spherical colloid subject to AC fields, accounting in full for the hydrodynamic and 
electrostatic interactions. The remainder of this article is organized as follows: In section 
[21 we give a brief introduction of our simulation model. We present the simulation results 
on the polarizability and effects of systematically varying the properties of the external 
AC fields in section [3l Finally, section [H concludes with a brief summary. 
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2. Model 

Our simulation system has three components: The solvent, the colloidal particle and 
the microions. The solvent is modeled as a fluid of DPD beads without conservative 
interactions, i.e., we only implement the DPD dissipative and stochastic interactions 
without conservative forces [S]. In the following, physical quantities will be reported in 
a model unit system of a (length) m (mass), e (energy), and e (elementary charge unit). 
In these units, the temperature of the system is kBT = l.Oe, the number density of 
the fluid is 3.0cr~^, and each solvent bead has a mass Im. The DPD friction coefficient 
is set to 'jDPD = 5.0y/me/a and the cutoff radius is r^ = 1.0a. The shear viscosity 
of the fluid was evaluated by analyzing pure Poiseuille flows in a microchannel by a 
method described in Ref. [27], and was found to be rjg = 1.23 ± O.Oly/me/a^, in good 
agreement with Ref. [25]. The diffusion constant of single solvent beads is measured to 
be Ds = 0.66 ± 0.02 a^/ej^. 

The colloidal particle is represented by a large sphere which has repulsive Weeks- 
Chandler- Anderson (WCA) [28j type conservative interactions with the fluid particles. 
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where the cutoff radius is set at the potential minimum re = r^ + 2^^^a. The radius 
of the colloid is R = rQ + a = 3.0a. To implement no-slip boundary conditions at the 
surface, a set of Ng DPD interaction sites is distributed evenly on the surface, with 
positions that are flxed with respect to the colloid center. These sites interact with the 
solvent beads through the DPD dissipative and stochastic interactions, using the same 
DPD friction constant as in the fluid, but twice the cutoff range. (The DPD cutoff has 
to be larger than the range of the WCA interaction in order to ensure that a sufficient 
number of fluid beads will interact with the DPD surface sites.) The total force exerted 
on the colloid is given by the sum over all DPD interactions with surface sites, plus the 
conservative excluded volume interaction, 

Ns 

Fc = Y.^i'\r.) + F^'^^^\ (4) 

1=1 
Here r^ denotes the position of z-th surface sites. Similarly, the torque exerted on the 
colloid can be written as 

Ns 

Tc = 5^Ff)(r,)x(r,-r_), (5) 

i=l 

where r^m is the position vector of the colloid's center-of-mass. Here the excluded volume 
interaction does not contribute since it points towards the colloid center. The force and 
the torque are used to update the position and velocity of the colloid in a time step 
using the Velocity- Verlet algorithm. The mass of the colloidal particle is M = lOOm and 
the moment of inertia is / = 360mcr^, corresponding to a uniformly distributed mass. 
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In addition, a total charge Q may be assigned to the center of the colloidal particle. 
To keep the whole system charge-neutral, a corresponding amount of counterions is 
then added to the solution. Salt microions may also be added as pairs of monovalent 
positively and negatively charged beads. We only consider the monovalent case where 
counterions or salt ions carry a single elementary charge ±le. Charges interact by 
Coulomb interactions, and the Bjerrum length Is = e^ / {AnemkBT) of the fluid is set 
to l.Oo". Apart from the electrostatic interaction with other charged particles, the ions 
have a short-range repulsive WCA interaction with the solvent beads and each other, 
which has the form of Eq. ([3]) with tq = 0. The conservative interaction with the colloid 
is also given by Eq. ([3]) with tq + o" = S.Oo". 

Specifically, we consider neutral colloids {Q = 0) and colloids with positive charge 
Q = +50e in salt-free solutions and in salt solutions with salt ion number density 
Ps = O.lo""^. In these systems, the diffusion constant of single microions is smaller than 
that of solvent particles and concentration dependent, i.e., Dj = 0.59 ± 0.02 a^/eJm in 
systems containing only counterions, and Dj = 0.47±0.01 a^fejm in systems containing 
counterions and salt. Figure [1] shows a representative snapshot of a single charged 
colloidal particle in salt-free solution. 




Figure 1. Snapshot of a colloidal particle in a salt-free solution. The surface sites are 
represented by the blue beads, and the white beads are counterions. Solvent beads are 
not shown here. 

Our simulations were carried out using the open source package ESPResSo [29], with 
a slight modification that allows us to incorporate a time-dependent external electric 
field. A cubic simulation box of linear size L = 20a with periodic boundary conditions 
was used for all simulations. Electrostatic interactions were calculated using the Particle- 
Particle-Particle Mesh (P3M) method [30], and a time step of At = O.Ola^ym/e was 
used for the integration. 
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As a test of our colloidal model, we have performed simulations of an uncharged 
colloid in a salt-free solution and measured the autocorrelation functions. Two 
autocorrelation functions were obtained from the simulations: the translational and 
rotational velocity autocorrelation functions 

(v(0) ■ v(t)) 



a(t) 



Mo)-^(^)) 



(6) 
(7) 



where v(t) and ci;(t) are the translational velocity and rotational velocity for the particle 
at time t, respectively. Figure [2] shows the simulation results averaged over ten runs 
with different random number generator initializations. 
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Figure 2. Translational (top) and rotational (bottom) velocity autocorrelation 
functions. The measurement is done for a single uncharged colloid with radius R = S.Oct 
in a salt-free solution. The temperature is ksT — l.Oe. 



For short time lags, both autocorrelation functions show exponential relaxation. 
The decay rate can be calculated using the Enskog dense-gas kinetic theory [3ll [32] . 
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They are plotted as blue lines in Figure [2] and show reasonable agreement with the 
simulation data when t < O.lcryrn/^. Similar results have been obtained using MPCD 
simulations [33 j . 

For long time lags, hydrodynamic effects set in and lead to a slow relaxation for 
autocorrelation functions of the colloidal particle [3l]. This so-called long-time tail is the 
manifestation of momentum conservation, as the momentum must be transported away 
from the colloids in a diffusive manner. Mode-coupling theory predicts an algebraic 
f-3/2 bgjiavior at long times for the translational velocity and t^^/^ for the rotational 
velocity [35] • The simulation results are plotted as green lines in Figure [2l The data 
are consistent with the theoretical prediction for t > lOay^m/e, but the rotational 
autocorrelation function exhibits large fluctuations for large times. This is mainly due 
to the fact that the statistics for long time values becomes very bad, and very long 
simulations are required in order to obtain accurate values. 

From the velocity autocorrelation function, the diffusion constant of the colloids 
can be calculated using the Green-Kubo relation 

D = l^°°rft(v(0)-v(t)). (8) 



In our case, this gives the diffusion constant D = 0.013 ± 0.002a ^^/e/m, which is in 
good agreement with the theoretical value for the diffusion constant of a Stokes sphere, 
D = ksT / {<6Txr]sR) = 0.0144 ± 0.0002a V'^T^. 

3. Results and Discussion 

In this section, we report our results for single colloidal particles of radius R = 3.0a 
under alternating electric fields. 

3.1. Displacement and velocity of the colloid 

We start with a system containing one colloidal particle of charge Q = +50e, an equal 
amount of counterions with charge — e, and monovalent salt ions with a number density 
O.lcr"^, resulting in a total of 450 negative and 400 positive microions in the solution. 

Figure |3] shows the colloidal response to an AC field in x-direction with frequency 
/ = 0.01^/e/m/a and amplitude Eq = 1.0e/{ae). After the AC field is turned on, the 
system is allowed to evolve for 10^ time steps. Then the measurement is done every 10 
time steps over a period of 10® time steps. In Figure [3l the time evolution of the electric 
field (top panel. Figure [3]^a)) is contrasted with that of various quantities of interest 
(Figure [3] (b)-(e)). We will begin with discussing the displacement and velocity, shown 
in Figure [3] (b) and (c). 

The electric field strength in this particular simulation was chosen relatively large, 
such that the variations of the position and velocity are above the thermal noise. The 
period of the applied field is T = 1// = 100a^^/rn/e, which is less than the diffusion time 
scale of the colloidal particle (r]^ = R^ / D ^ QAOa-\/m/e, where D = 0.013a ^/ejm 
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Figure 3. Dynamics of a ciiarged colloidal particle with charge Q = +50e and radius 
R = 3.0(7 in a salt solution (number density of salt ions 0.1a~^) under the influence 
of an AC field with frequency / = O.Ol^/eJm/a and amplitude Eq = l.Oe/{ae). 
(a) Instantaneous electric field, (b) Position of the particle along the field direction, 
(c) Translational velocity along the field direction, (d) Effective charge by counting 
number of ions inside the radius r < 4. Oct. (e) Dipole moment in the field direction. 
The effective charge and the dipole moment are discussed in section [321 



is the diffusion constant of the colloidal particle). On the one hand, the particle thus 
experiences many oscillation before diffusing over the distance of its radius. On the 
other hand, the period is still much larger than the diffusion time scale of the fluid 
(r)) = a'^/Ds ~ l-'bo^frnje^ where Ds = 0.66a^^/m/e is the diffusion constant of 
solvent particles). 

To rationalize the results, one may thus adopt a simplified picture where the 
colloidal particle is viewed as a charged sphere immersed in a viscous fluid in an AC-field. 
The equation of motion can be written as 

Mx + -fx- QEoe''^' = 0. (9) 

The friction term is inversely proportional to the velocity and its value can be 
approximated using the Stokes friction 7 = GnrisR ^ TOy/me/a. We solve the differential 
equation ([H]) by assuming that the position takes the form x = Ae*'^*, where A is a 
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complex number that characterizes the magnitude and phase of the oscillation. The 
results are 

\A\ = , ^^° , Arg(A) = tan-^(-^)--. (10) 

Inserting the parameters of our system, we obtain that the amplitude of the position 
oscillation should be about llcr, and the phase should be advanced by about 7r/2 
compared to the applied field. Also, the amplitude of the velocity oscillation should 
be O.TlA/i/Vri and the velocity should be almost in phase with the applied field. 

The simulation results shown in Figure [3] are in good qualitative agreement with the 
simplified theory. The position variation clearly exhibits a 7r/2 phase shift ahead of the 
applied field, while the velocity profile is in phase. However, the theory overestimates 
the amplitude of the oscillation. The actual amplitude of the position oscillation is about 
5cr, and that of the velocity is about 0.3ye/m. One reason for this discrepancy is that 
counterions condense on the colloids surface, which effectively reduces the total charge 
of the colloids. The binding of the counterions and the colloid also increases the effective 
mass of the complex, which further reduces the amplitude of the oscillations. Another 
effect comes from the no-slip boundary condition of the colloidal particle: The fluid 
particles adjacent to the colloid surface are dragged along with the colloid, effectively 
increasing its mass. 

One important quantity that characterizes the colloidal response to AC-fields is the 
mobility /i(cj). In the presence of an external field Ee*"^*, the particle oscillates with the 
velocity Ue"^* with 

U = /iME, (11) 

where the mobility /i(u;) is a complex function of the frequency. In the simulation, the 
time sequence of the velocity is noisy due to the thermal fluctuation. To obtain a clear 
signal and extract the mobility, we apply Fourier transform to the time sequence of the 
velocity. Figure S] shows the Fourier transform of the time sequence in Figure [3t^c). The 
real and imaginary parts correspond to the in-phase and out-of-phase components of 
the signal with respect to the applied field. They both exhibit peaks at the frequency 
of the external electric field, and the peak values divided by the electric field strength 
give the real and imaginary part of the mobility at that frequency. 

3.2. Charge distribution and dipole moment 

Apart from allowing one to study the displacement and velocity of the colloidal 
particle, the simulations also provide valuable information about the electric double 
layer surrounding the colloid. During the simulation, we have collected the histograms 
of counterion and colon distances to the center of the colloidal particle. Figure [5]^a) 
shows the time-averaged charge distribution around the colloidal particle. Due to the 
strong Coulomb interaction with the positively charged colloids, counterions accumulate 
close to the colloidal surface while the colons are depleted in the same region. The 
counterion concentration is peaked around r = 3.2a from the colloid center. Beyond 
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Figure 4. Fourier transform of tire time sequence of the colloid velocity of Figure [3]Jc). 
The top panel shows the real part corresponding to the in-phase component of the 
velocity, and the bottom panel is the imaginary part corresponding to the out-of-phase 
component. Both the real and the imaginary parts feature peaks at / = 0.01 ye/m/cr, 
the frequency of the external field. 



that distance, the counterion concentration decays rapidly to the bulk value, which is 
reached at around r = 6. Oct. At the same time, colons are barely present within the 
distance r < A.Oa, and the concentration slowly increases to the bulk value, which is 
also reached at around r = 6.O0". 

The ion cloud surrounding a charged colloidal particle is commonly grouped into 
two regions separated by a "plane of shear", which is defined as a hypothetical plane 
that separates an immobilized liquid layer and the mobile fiuid. In the region inside the 
shear plane, also called the stagnant layer, the ions are strongly bound to the surface and 
have limited mobility. Due to the friction between the liquid and the particle surface, 
ions inside the stagnant layer follow the surface. Outside the stagnant layer, the ions are 
relatively free to move. The outer region up to the point where the ion concentration 
reaches the bulk level is called diffuse layer. The thickness of the diffuse layer is roughly 
determined by the Debye screening length, although it may effectively be smaller due to 
nonlinear effects already at moderate surface charges [36]. The Debye screening length 
is given by 
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(12) 



where Zi and nj(oo) are i-th ion's valence and bulk number density, respectively, and 
the summation is over all ion species. In our system, the Debye screening length is given 
by Id = 0.89cr in pure solution (neglecting the small contribution of the counterions) . 

Determining the extension of the stagnant layer is important because it contains 
the ions that contribute to the effective charge of the particle. In simulations, however, 
it is difficult to pin down the exact position of the shear plane. Different methods have 
been proposed in the literature [371 EH ISH HQl ST] . One popular method is the inflection 
criterion [371 EHl |39]. Here one flrst computes the integrated counterion distribution 
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Figure 5. (a) Average number density of counterions and colons as a function of the 
distance to the colloid center. The average has been taken over 10^ time steps, (b) 
The integrated counterion distribution /(r) used to determine the cutoff radius for the 
effective charge. The inflection point is identified as the zero of the second derivative 
of I{r) with respect to 1/r. 



I{r) by counting the number of counterions inside a radius r with respect to the colloid 
center, 

I{t) = ^ j\^{T)dV, (13) 

where p-{r) is the counterion number density and the normalization constant A^_ is the 
total number of counterions. The value of /(r) represents the fraction of counterions 
which are within a distance r from the colloid center. The cutoff radius for summing 
up the effective charges is then set at the inflection point (the point where the second 
derivative is zero) of the integrated counterion distribution plotted versus the logarithmic 
distance to the center. Although this criterion has originally been devised for the 
static case, here we will attempt to use it for the dynamic situation by applying it 
to the averaged counterion distribution. Figure Mjo) shows the integrated counterion 
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distribution. The inflection point is identified at r = 3.9cr. 

Another method to determine the shear plane is to compute the colloid-ion velocity 
correlation HOl ST] 



Cci{r,T) = {vc{t)vi{r,t + T)) - {v,{t)) (t;^ (r, t + r) ) , (14) 

where Vc{t) is the velocity of the colloidal particle and f j(r, t) is the velocity of ions which 
is located at r with respect to the colloid center. Since the motion of the colloid and the 
ions are strongly perturbed in the direction of the external field, we only consider the 
correlation in the other two directions. Figure [6] shows the correlation functions for five 
different r. For ions which are located in the nearest surroundings of the colloid surface, 
the correlation in the limit of r — )■ is positive, indicating that those ions are moving 
along with the colloid. At larger distance, the ion motion becomes anti-correlated with 
that of the colloid, and the correlation Cci{T — )■ 0) turns negative. The position where 
the correlation first changes sign can also be viewed as a definition of the shear plane. 
Using this method, one locates the shear plane at r = 3.8cr, which is close to the value 
determined with the infiection criterion. 
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Figure 6. Colloid-ion velocity correlation function Cci{r,T). The values are averaged 
over two directions normal to the external field. 



In order to characterize the dielectric response of the colloidal particle under 
external fields, we must calculate the dipole moment of the system. Whereas only 
the ions inside the stagnant layer contribute to the effective charge, the dipole moment 
also has significant contributions from the diffuse layer. Figure O^a) indicates that the 
ion concentrations reach their bulk values at approximately r = Q.Oa, therefore we 
computed the dipole moment using a cutoff of Q.Oa from the colloid center. The time 
sequence of the effective charge and the dipole moment are plotted in Figure M^d) and 
(e), respectively. The averaged effective charge is roughly +15e, corresponding to about 
condensed 35 counterions on the colloidal surface. The dipole moment oscillates with 
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the same frequency as the apphed field. The complex polarizability a{uj) is obtained by 
applying a Fourier transform to the time sequence of the dipole moment, similar to the 
complex mobility discussed in last section. 

3.3. Influence of the strength of AC fields 

Most theories on the dielectric response of charged colloids are based on the linear 
response theory and assume the external electric field to be weak. In our simulations, 
we can apply large electric fields to examine the effect of the field strength and possible 
nonlinear effects. In this subsection, the frequency is set at / = O.Ol^/e/m/a, which is 
in the low-frequency regime (see next section). The field strength is varied from 0.01 to 
lOe/i^ae). By covering a wide range of electric field amplitudes, we are able to estimate 
the limit of validity of the linear response theory. Figure [7] shows the amplitudes of the 
colloid velocity and dipole moment as a function of the field strength. 
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Figure 7. (a) The amplitude of tlie colloid velocity as a function of the field strength, 
(b) The amplitude of the dipole moment as a function of the field strength. The 
frequency of the applied field is set at / = 0.01^ye/m/a. 



In weak fields, the velocity grows linearly as a function of the field strength up 
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to Eq ^ 2.0e/ (ere), corresponding to a constant mobility. For higher field strength, the 
growth slows down, which implies that the mobility decreases. This is different from the 
case of constant electric field, where the mobility is found to increase above the linear 
region [H]. Alternating electric fields of large field strength can excite colloidal motion 
with signatures of various frequencies besides the frequency of the applied field. This 
can be seen from the Fourier transform of the colloid velocity, shown in Figure [H where 
peaks at higher frequencies appear for high field strength Eq = 10e/{ae) (see Figure H] 
for comparison). 
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Figure 8. Fourier transform of tire time sequence of the colloid velocity for a strong 
external field with frequency / — O.Qly^e/m/a and amplitude Eq = 10e/(iTe). 



Looking at the dipole moment, one finds that the deviation from the linear regime 
sets in at about Eq ^ 1.0e/{ae). Compared to the colloid velocity, the deviation from 
the linear dependence is much more pronounced for the dipole moment. This is due to 
the fact that the dipole moment measures how ions are distributed around the colloidal 
particle, and the electric double layer can easily be infiuenced by the external field. 
For intermediate field strength, the main effect of the external field is to deform the 
shape of the electrical double layer, as the ions in the diffuse layer are relatively mobile. 
The shape of the ion clouds is elongated in the field direction. Figure |9] shows the ion 
distribution for three different field strength, Eq = 0.1, 1 and 10£:/(o"e). In the weak 
field region {Eq = 0.1 to 1.0e/{ae)), the change of ion distribution is small. When 
the field strength increases to 10e/{ae), the deviation becomes notable. The peak of 
the counterion distribution is reduced, and a second peak of colons close to the colloid 
surface develops. 

When the field is strong enough, it may disturb the ions inside the stagnant layer. 
Some of the ions originally bound to the colloidal surface will be ripped off at large 
enough electric field. This is demonstrated in Figure M where the counterion peak 
is reduced at higher field strength. The required field strength can be estimated by 
equating the force exerted by the external field and the electrostatic force from the 
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Figure 9. Comparison of tlie ion distribution for three different field strength 
Eo ~ 0.1,1 and 10£/(fTe). The top and bottom figures are for the counterion and 
coion distribution, respectively. 



charged colloids. 

Eq 



Q 



eff 



AnemR'^e 



1.7e/{ae), 



(15) 



where we have used Qes = 15e. This number seems to determine the point where the 
dipole moment as a function of field strength starts to level off, indicating the breakdown 
of the weak field assumption. 



3.4- Frequency dependence 

Finally, we investigate the effect of varying the frequency of the external AC field. The 
amplitude of the field is chosen in the linear region, Eq = 0.5£:/(cre). The frequency is 
varied from / = 10~^ to l.^^elmjo. A separate simulation with a constant electric field 
is also performed to provide the reference at low- frequency limit. As in most previous 
simulations, the solution has a salt ion number density of O.lcr"^, corresponding to 
the Debye length Id = O.S9a in a system without counterions and Id = O.S7a in a 
system with 50 counterions. Thus the Debye length is smaller than, but comparable 
to the colloidal radius {R = 3. Oct). We compare the simulations with the prediction of 



the Maxwell- Wagner-O'Konski (MWO) theory, which is briefly sketched in Appendix A 



Even though this theory was originally developed for situations where the electric double 
layer is much thinner than the colloid, we shall see that it still captures the most 
important features of the frequency dependence on a semi-quantitative level. 

We start with the easiest case of an uncharged colloidal particle, which has no 
electric double layer. Figure [10] shows the polarizability for an uncharged colloidal 
particle in salt solution with ions number density of O.lo""^. The real part of the 
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polarizability shows the in-phase component of the dipole moment with respect to 
the external field, while the imaginary part gives the out-of-phase contribution. The 
polarizability in the low-frequency limit has a negative real part, indicating that the 
induced dipole is anti-parallel to the apphed field. Uncharged colloids are not surrounded 
by an electric double layer, and the dipole moment is induced solely by the motion of 
the salt ions. The positive ions move in the direction of the electric field, and tend to 
accumulate on the back side of the colloid, while the negative ions accumulate on the 
front end of the colloid. This leads to an effective dipole moment which points in the 
opposite direction than the external field, sometimes referred to as volume polarization 
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Figure 10. Real and imaginary part of tire complex polarizability of an uncharged 
particle as a function of the frequency of applied electric field. The field strength is set 
in the linear region _Eo = 0.5£/(ae). The number density of the salt ions in the solution 
is O.Ict""^. The points are simulation results. The solid lines give the prediction from 
the Maxwell- Wagner theory with the effective colloid radius i?eff = 2.93cr. The dashed 
line shows revised prediction by taking into account of inertia effects. See text for 
explanation. 



The solid line in Figure [TO] shows the prediction from the Maxwell-Wagner theory. 
The limiting behavior of the polarizability at frequency / — ;■ was obtained by an 
independent run in a DC field and used to calibrate the effective colloidal radius R via 
a(0) = AnemKoR^, which reduces to a(0) = —R^/2 in our system (see Appendix A ). 
We obtain R^g = 2.93 ± 0.07cr, which is reasonably close to the "physical" radius 
R = 3.0a. The main sources of discrepancy are the soft character of the excluded 
volume interactions and the small size of the colloids. Having set the effective colloidal 
radius, we can calculate the Maxwell- Wagner prediction without further fit parameter, 
and the resulting curve is in quite good agreement with the simulation data (Figure 
[TOl solid line). Most notably, the theory predicts a crossover between two regimes 
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which is recovered in the simulation at roughly the predicted crossover frequency. The 
simulations feature a slight overshoot of the polarizability in the crossover region, which 
can be explained by considering the inertia effect of the microions, which are not entirely 
negligible in our simulations. (In real systems, they are negligible.) When microions 
with mass are immersed in a viscous fluid, their response to the external electric fields 
depends on the frequency. Thus the conductivity becomes frequency-dependent. The 



inertia effect can be incorporated into the theory (see Appendix A) and the theory can 
be revised accordingly, giving the dashed line in Figure [THl It predicts small overshoots 
in both the real and the imaginary part of the polarizability, which are consistent with 
the simulation. 

The situation is more complicated when the colloids are charged, due to the presence 
of the electric double layer. Figure [TT] shows the polarizability for a charged colloid 
[Q = 50e) in salt solution (same salt concentration as before). 
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Figure 11. Same as Figure [10] for a charged colloid {Q — 50e). The points are 
simulation results. The solid lines give the prediction from the MWO theory, using an 
effective colloid size Ras = 2.93 cr and the surface conductance Ka — O.^Ae'^/ (ay/me). 
The dashed line shows revised prediction by taking into account of ions' inertia effect. 
See text for explanation. 



In the low- frequency region / < IO^^a/^/^/o", the external perturbation is slow 
enough that the system can follow. The dipole moment is in phase with the applied 
field, as reflected by the fact that Im{a} is close to zero. The real part Re{a} is positive, 
in contrast to that obtained for uncharged colloids, indicating that the dipole moment 
points in the same direction as the applied field. In the absence of external fields, the 
center of the counterion cloud and the center of the spherical particle coincide, and the 
dipole moment is zero. If one applies an electric field, the positively charged particle 
moves in the direction of the field, while the negative ion cloud moves in the opposite 
direction. The resulting net dipole moment points in the same direction as the external 
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field, therefore one lias Re{a} > in the low-frequency limit. In the opposite limit of 
high frequency, the colloid and the ion cloud can no longer follow the field, thus both 
Re{a} and Im{a} converge to zero. At intermediate frequencies / ~ 0.1(e/mY^'^a~^, 
the real part Re{a} crosses over from positive to zero with an overshoot below the 
transition frequency and a slight undershoot to negative values after the transition 
frequency. The imaginary part Im{a} assumes a maximum, indicating that the response 
is out of phase and that there is high dissipation. 

According to the Maxwell-Wagner-O'Konski (MWO) theory, the effect of the 
electric double layer can be taken into account by introducing a surface conductance 
in the Maxwell- Wagner theory. Here we have determined the value of the surface 
conductance by calculating the / — )■ limit of the polarizability from an independent 
simulation in a DC field, assuming that the effective radius is the same for charged 
and uncharged colloids. The simulation gives the surface conductance K^ = 0.24 ± 
0.07e'^/{a^/me), which is consistent within the error with the theoretical estimate based 
on the extended MWO theory [9] (see appendix), K„ = 0.19e^/((T\/me). Using this 
value, we can again calculate the theoretical polarizability as a function of frequency 
without any fit parameter. The resulting curves without and with ion inertia effects 
are shown in Figure [11] with solid and dashed lines, respectively. The prediction of 
the MWO theory is in surprisingly good agreement with the simulations, given how 
questionable its validity is. The theory captures the main qualitative features, and 
roughly the correct crossover frequency. Similar to the uncharged case, the dielectric 
response at high frequency are affected by the finite mass of microions. The revised 
theory shows better agreement with the simulation after the transition frequency. The 
overshoot in the real part Re{Q;} may be an indication of the low-frequency dielectric 
dispersion [10], and further studies are required to clarify this problem. 



4. Summary 

We have carried out mesoscopic molecular dynamics simulations for a colloidal particle 
under alternating electric fields. We have accounted in full for the hydrodynamic and 
electrostatic interactions as well as thermal fluctuations, using Dissipative Particle 
Dynamics and a Particle-Particle-Particle Mesh method for evaluating Coulomb 
interactions. 

We obtained information about the displacement and velocity of the colloidal 
particle. A simple classical theory is proposed to explain the simulation results at 
low frequency and shows positive agreement. We have also considered the contribution 
from the electric double layer surrounding the colloid. Two different methods, namely, 
the inflection criterion and the colloid-ion correlation function, were used to determine 
the position of the shear plane, and it was shown that they both give the same results. 
Furthermore, the dipole moment was calculated by integrating all contributions of ions 
from the colloid surface to the point where the ion concentration reaches the bulk value. 

The effect of the external field strength was investigated systematically. By 
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computing the mobility and the polarizabihty for various field strength, we have 
examined the validity of the linear response theory. Nonlinear effects were found to 
set in for electric strength up to 1.0e/{ae), which is higher than the values reported for 
constant electric fields [HI [191 US]- This can be contributed to the excitation of high- 
frequency modes and the stripping of the counterions that are bounded to the colloid 
surface. 

Furthermore, the frequency dependence of the polarizabihty was studied for an 
uncharged colloid and a charged colloid in a salt solution. The simulation results were 
compared with the predictions of the Maxwell- Wagner-O'Konski theory. Despite the fact 
that thermal fiuctuations are substantial and that the size of the colloids is comparable 
to the Debye screening length, the MWO theory performed surprisingly well, especially 
for uncharged colloids. The inertia effect of microions can also be incorporated into the 
theory to explain the high frequency behavior. For charged colloids, the theory and 
simulations show noticeable deviations, indicating that the thin assumption of a thin 
electrostatic double layer is questionable for our system. 
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Appendix A. Maxwell- Wagner-O'Konski theory 

In this appendix, we give a short introduction to the Maxwell- Wagner-O'Konski theory. 
More detailed information can be found in Refs. [6l [71 El |9] . 

When loss is present, the dipole moment of a particle immersed in fluids exhibits 
a phase lag with respect to the external sinusoidal fleld. A complex effective dipole 
moment can be written as p = A7remK{u)R^ E (cf. Equation ([2|)), where the Clausius- 
Mossotti factor K{uj) is a complex number containing both the magnitude and the phase 
information about the effective dipole moment. In the Maxwell- Wagner theory, it has 
the form 

^«>4) = |^, (A.l) 

where e* and ej^ are the complex dielectric constants of the particle and the medium, 
respectively. They are deflned as 

K K 

^p — ^p ~r -. — ; ^m — ^m-\ : ) l^"^J 
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where e (without the star) and K are the permittivity and conductivity, respectively. 
For e.^ = Cp = e, the Clausius-Mossotti factor can be rewritten as 

K-k/j^ (A.3) 

with 

_ -fvp ~ J^m 1 _ _ o^ ( A A"] 

In our system, we have Cp = e^ = l/47r, i^p = 0, thus Kq = —1/2 is negative. The 

conductivity of the medium K^ is related to the microion diffusion constant Dj via 

Km = nio^e^Di/{kBT) for 1-1 electrolytes. 

The classical Maxwell- Wagner theory fails to explain the dielectrophoretic 

properties of latex particles. Latex has a low intrinsic conductivity, but the measurement 

indicated that the particle conductivity is high. It was concluded that the electric 

double layer surface significantly contributes to the particle conductivity. This was first 

demonstrated by O'Konski |8j 

2K 
Kp^Kp + ^ (A.5) 

ri 

where K„ is the surface conductance (unit S instead of S-m~^ for conductivity) due 

to the electrostatic double layer. The extra contribution to the conductivity increases 

the Maxwell- Wagner relaxation frequency without introducing an extra dispersion (no 

extra extrema in the imaginary part of the polarizability). The idea is to encapsulate 

the contribution from the ion cloud around the colloids into one single parameter (the 

surface conductance). This is justified for high ionic strength, highly charged colloids 

and high frequencies. Another important effect of K^j is to shift the prefactor Kq in Eq. 

f lA.3|) towards positive values, to the point that the dielectric response changes sign if 



K^ is sufficiently large. 

We still need a way to obtain the surface conductance K^j from the effective surface 
charge density a in the simulation. This can be done by the procedure outlined in [9]. 
First, the ^-potential is solved using the following formula [TT] . 



(^mkBT 

a - 



el 



D 



2 sinh ( — ; — — ■ I -|- ^ ,, tanh 



2kBTj R/Id \4:kBT 



(A.6) 



The surface conductance then can be related to the ^-potential by Bikerman's expression 



K, = Id 



'^^Pflal^l) - ' 



(1 + 3m)Km, (A.7) 



where m is a dimensionless ionic drag coefficient [HI US]. Bikerman's expression is 
applicable for 1-1 electrolytes, and the contribution from colons has been neglected. 
The ionic drag coefficient is related to the viscosity rjs and ion diffusion constant Dj, 
2em{kBTf 



m 



'irjse^Di 

{m = 0.092 for our system, and m 
obtain the surface conductance K„ 







(A.8) 


0.18 for KCl). With Q = 


= 50e and Rcs 


= 2.93(T, we 


= 0.19eV(o^v/^). 
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The inertia effect of ions can be taken into account by considering a charged particle 
with charge e, mass m immersed in a viscous fluid. The equation of motion for the 
particle under an AC field E^e^'^^ is 

mx = -^x + eEoe''^\ (A. 9) 

where 7 = ksT/Dj is the friction constant. After solving the equation of motion, one 
finds that the velocity of the particle has the form 

l-z(^) 

Since the conductivity of salt solutions is proportional to the ion's velocity, it is 
reasonable to assume that the conductivity has the same frequency dependency 



The dashed hues in Figure [TO] and [TT] are obtained by substituting (lA.lip into the 
Clausius-Mossotti factor ( 1A.3I) . 
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